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Abstract 

We present analytical potential-density pairs in three dimensions 
for the gravitational field of galaxies, obtained by thickening the multi- 
polar expansion up to the quadrupole term. These may be interpreted 
as generalizations of the Miyamoto and Nagai potential-density pairs. 
With a suitable restriction on the possible values of the multipole mo- 
ments, the density distributions are positive and monotone decreasing 
functions of the radial and axial coordinates. 
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1 Introduction 

There are several three-dimensional analytical models in the literature for 
the gravitational field of different types of galaxies and galactic compo- 
nents. Jaffe (1983) and Hernquist (1990) discuss models for spherical galax- 
ies and bulges. Three-dimensional models for flat galaxies were obtained 
by Miyamoto and Nagai (1975) and Satoh (1980). de Zeeuw and Pfen- 
niger (1989) considered a set of ellipsoidal models appropriate to galactic 
bars. Long and Murali (1992) derived simple potential-density pairs for a 
prolate and a triaxial bar by softening a thin needle with a spherical and 
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a Miyamoto and Nagai potential, respectively. See Binney and Tremaine 
(1987) for a discussion of other galactic models. There also exist several 
general relativistic models of disks, e.g., Morgan and Morgan (1969, 1970), 
Bicak et al. (1993), Lemos and Letelier (1994), Gonzalez and Letelier (2000, 
2004), Vogt and Letelier (2003, 2005a). Recently, a general relativistic ver- 
sion of the Miyamoto and Nagai models was studied (Vogt and Letelier 
2005b). 

The potential-density pairs obtained by Miyamoto and Nagai (1975) 
are inflated versions of the thin-disk family first derived by Toomre (1963). 
In this work we consider a set of three-dimensional potential-density pairs 
obtained by using the same procedure of Miyamoto and Nagai on the mul- 
tipolar expansion up to the quadrupole term. This generates a sequence 
of potential-density pairs that reduces to that of Miyamoto and Nagai for 
particular values of the multipole moments. This is done in subsection 12.11 
and subsection 12.21 The thin-disk limit is investigated in subsection 12.31 
and corresponds to generalizations of Toomre's family of disks. Finally, the 
results are discussed in section |3 



2 Multipolar Models for Flattened Galaxies 

A general expression for a multipolar expansion in spherical coordinates can 
be cast as 

n=0 

where P n are the Legendre polynomials and a n are coefficients related to the 
multipolar moments. We consider only the expansion up to the quadrupole 
term (n = 2). In cylindrical coordinates the explicit form reads 

Gm Dz Q(-R 2 + 2z 2 ) 

VWT^ {R 2 + z 2 fl 2 2{R 2 + z 2 fl 2 ' ( ' 

where D and Q are the dipole and quadrupole moments, respectively. 

In the following three-dimensional models the mass-density distribution 
is obtained directly from Poisson equation, 

' r *,RR + ^r + *,zz). (3) 



r 4ttG V ' R 

Other physical quantities of interest are the circular velocity, v c of particles 
in the galactic plane, the epicyclic frequency, k, and the vertical frequency, 
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z/, of small oscillations about the equilibrium circular orbit in the galactic 
plane. They are calculated with the expressions (Binney and Tremaine 
1987): 

vl = R*R, (4) 
3 

= $,RR + ^,R, (5) 

v 2 = (6) 

where all quantities are evaluated on z = 0. The stability conditions are set 
by k 2 > and v 2 > 0. 

2.1 Generalized Miyamoto and Nagai Model 2 



We first consider Q = 0, and apply the transformation z — ► a + \/z 2 + & 2 
on the multipolar expansion equation P|). where a, 6 are non-negative con- 
stants. Using equation ©, we obtain 



i? 2 + (i + e) 2 



7/2 



{i? 4 (i - b) + i? 2 (i + o [d(i - so 



+(1 + 0(2 + 30] + (1 + 3 [2^(1 + 40 + (1 + 0(1 + 30] } , (7) 

where the variables and parameters were rescaled in terms of a: R = R/a, 
z = z/a, b = b/a, D = D/(Gma), p = mp/a 3 , and £ = \j z 2 + b 2 . For 
the particular value D = 1, equation Q reduces to the mass density of the 
Miyamoto-Nagai model 2. 

Unfortunately the density distribution, equation ( |7| l. has some defects. 
For certain ranges of the parameters b and -D it is not a monotone decreasing 
function of the radial and axial coordinates, and even has domains with 
negative densities. To overcome this, we impose a reasonable restriction that 
the derivative of the density distribution with respect to the R coordinate 
along z = should have only one critical point at (i?, z) = (0, 0) as well 
as the derivative with respect to z along R = 0. A graphical study of the 
resulting equations is displayed in figures ^i-b. In figure^, the curves of 
p £ = are plotted as functions of the radial coordinate and the parameter 

D for some values of b. In this case decreasing b narrows the interval of D 
for which there are no other critical points, as on R = 0. Figure ^ shows 
the curve of p^ = as functions of z and D for the same values of b. One 
can see from both graphs that the allowed interval for the parameter D for 
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a fixed value of b in the case (a) is contained in the allowed interval in case 
(b). 

The expressions for the circular velocity, epicyclic frequency and vertical 
frequency follow from equations Q- © 



R 2 



R 2 + (l + b)(l + b + 3D) 



R 2 + (l + b) 2 
1 



R 2 + (1 + If 



7/2 



5/2 ' 

| J R 4 + J R 2 (l + 6) +/))- :\D 



(8) 



-4(1 + 6) 3 (1 + 6 + 31))}, 



, 2 _ R 2 (l + b-D) + (l + b) 2 (l + b + 2D) 



R 2 + {1 + b) 2 



5/2 



(9) 
(10) 



where v 2 = Gmv 2 /a, k 2 = GmK 2 /a 3 , and v 2 = Gmv 2 /a 3 '. The condition 
v 2 > and the stability conditions k 2 > and ^ 2 > also impose restrictions 
on the possible values of D and b. We find that D > — (1 + 6)/3 ensures the 
positivity of the circular velocity and of the square of the epicyclic frequency, 
whereas the square of the vertical frequency is always non-negative if D > 
-(l + 6)/2. 

In figures EK~c we plot some isodensity curves of the density function, 
equation 0, with parameter b = 0.5 and (a) D = 1, (b) D = 0.5, and (c) 
D = —0.1. For these values the density is a monotone decreasing function, 
as can be checked from figure ^ It is clearly seen that as the parameter 
D decreases, the density distribution becomes more flattened. Figure EJi 
and figures ^-b show, respectively, curves of the velocity profile, equation 
(jHJ, the epicyclic frequency, equation ©, and of the vertical frequency, 
equation ()10j) . for the same parameters as in figures |^ _c - With decreasing 
D the radius where the highest circular velocity occurs is increased and the 
epicyclic frequency is lowered. The vertical frequency becomes lower near 
the center, but for R <; 2 it is increased. 
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Figure 1: (a) Curves of p ^ = as functions of R and D for b = 2, 1, 0.5, 
and 0.1. (b) p : z = as functions of z and D for b = 2, 1, 0.5, and 0.1. 
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Figure 2: Constant-density curves of equation Q with parameters b = 0.5 
and (a) D = 1, (b) D = 0.5, and (c) D = —0.1. (d) The circular velocity v c 
(equation (jHJ)), in the galactic plane for cases (a)-(c). 




Figure 3: Curves of (a) the epicyclic frequency, k (equation @), and (b) the 
vertical frequency, v (equation HlOjO . for the same parameters as in figure |21 
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2.2 Generalized Miyamoto and Nagai Model 3 

We now consider the full expression equation (J2J) and apply the transforma- 
tion z — > a + \/ z 2 + b 2 . The resulting mass density distribution reads 

b 2 



8tt£ 3 



i? 2 + (l + 2 



9/2 



|2fl 6 (l -Z>) + 3-R 4 2(1 + 3 - 6Z?£(l+f) 



-3Q] + 3i? 2 (l + f) 2 [2(1 + £) 2 (1 + 20 + 2D(1 + - Q(l + 25£) 
+2(1 + 4 [(1 + 2 (1 + 30 + 25(1 + 0(1 + 40 + 3Q(1 + 50] } , (11) 

where the variables and parameters have been rescaled, as in subsection 12. II 
and Q = Q/(Gma 2 ). For particular values D = 1, Q = 2/3 we recover the 
Miyamoto and Nagai model 3. 

Also here the range of the parameters must be restricted to produce 
physically acceptable density distributions. In figures EJi-b we show, re- 
spectively, some curves of p f, = along z = and p t s = along R = 

as functions of Q for some sets of values of D and b. The expressions for 
the circular velocity, epicyclic frequency and vertical frequency follow from 
equations 

_o R 2 



R 2 + (l + b) 2 



7/2 



|2i? 4 + R 2 



4(1 + by + 6D(1 + b) - 3Q 



-2(1 + bf [(1 + bf + 3D{\ + b) + 6Q] } 



(12) 



k 2 



i? 2 + (1 + 6) 2 



9/2 



{ 



2R 6 + 3i? 4 



4(1 + bf - 2D(1 + b) + Q 



+18R 2 (1 + bf [(1 + 6) 2 + D(l + 6) - 3Q 
\4 f/i , r\2 



*> 2 



-8(1 + 6f [(1 + 6) z + 3D(1 + 
1 



2b 



R 2 + (1 + 6) 2 



7/2 



6) + eg] } , 

{2i? 4 (l + b - D) + i? 2 (l + 6) [4(1 + bf 



(13) 



+2£»(1 + 6) - 9Ql + 2(1 + bf [(1 + bf + 2D(1 + b) + 3Q] } , (14) 

A graphical analysis of the constraints on the parameters imposed by v 2 > 0, 
k 2 > and v 2 > shows that they are contained in the restrictions imposed 
by the gradient of the mass density. 



7 



Q 





b=0.5, D=1.0 
6=0.5, 6=0.5 
6=0.5, 5=0 
b=0.5, 5= - 0.5 

6=0.5, 5= -1.0 
6=0.1, 6=1.0 
6=0.1, 6=0.5 
6=0.1, 6=o 
6=0.1, 5= -0.5 
6=0.1, 6=- i.o 



0.5 



-M- 



-0.5 




b=0.5, D=1.0 
6=0.5, 5=0.5 
6=0.5, 5=0.0 
.§=0.5, 6= -0.5 

6=«.5, 6= -i.o 
6=Ooi, 5=1.0 
6=o.t?5=o.5 
6=0.1, 6=o.o 

6=0.1, 6= %5 

fc=o.i, 6= -ft) 



Figure 4: (a) Curves of p ^ = as functions of R and Q for some values of 

the parameters b and D. (b) Curves of p^ = as functions of z and Q for 
some values of the parameters b and D. 
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Figure 5: Constant density curves of equation with parameters b = 0.5 
and (a) D = 1, Q = 2/3, (b) D = 1, Q = -0.1, (c) D = 0, Q = 2/3, and 
(d) D = 0, Q = -0.1. 
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Figure 6: Curves of (a) the circular velocity v c equation (|T2|) in the galac- 
tic plane, (b) the epicyclic frequency R equation (|T3|) and (c) the vertical 
frequency v equation (|14j) for the same parameters as in figure El 



Some level curves of the mass density, equation (fTTj) . are displayed in 
figures [5^,-d. In general, as D and Q are lowered for a fixed b, the mass 
distribution profile becomes flatter. We also note from figures HJi-c that the 
maximum of the circular velocity is shifted to larger radii, and the epicyclic 
frequencies and vertical frequencies are lowered near the center. 



2.3 Thin Disk Limit 

To investigate the thin limit of the potential-density pairs deduced in the 
previous sections, it is more convenient to rederive an expression for the 
surface density, a, by applying the transformation z — > a + \z\ on the mul- 
tipolar expansion, equation © and using the well-known relation (Binney 
and Tremaine 1987), 
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where the right-hand side is evaluated at z —* + . We obtain 
1 



a 



4vr(i? 2 + l) 7 / 2 



2(R 2 + 1)^(1 - D) + 6D(ir + 1) - 3Q(3i^ - 2) 



(16) 

The variables and parameters were rescaled as in subsection 12.21 and a = 
ma/a 2 . In particular cases with D = 1, Q = and D = 1, Q = 2/3 we 
obtain Toomre's models 2 and 3, respectively (Toomre 1963). Expressions 
for the circular velocity and epicyclic frequency are given by equation (|T2*|) 
and equation (|13|) with 6 = 0, and the vertical frequency is calculated from 
equation © evaluated at z — ► + : 



£ 2 



1 



2{R 2 + l) 9 / 2 L 



2R* + 9R 4 {2D -Q) 

+6R 2 (-1 + D + 12Q) - 4(1 + 3D + 6Q) 



(17) 



Figure [7| shows curves of <7 ^ = (lines without circles) and D = (lines 

with circles) as functions of Q and R for some values of the parameter D. 
In the thin limit, the disks always have regions of vertical instability. The 
Kuzmin disk (D = Q = 0) is stable for R > y/2; Toomre's model 2 (D = 1, 
0) is stable for ,R > v / 2\/-2 + \/6; if Q = the disk is stable for 



Q 



R > 



V2 



1-9D + J3(27D 2 + 2D + 3) 



1/2 



(18) 



if D > -1/3; and 



< R < 



V2 



1-9D- J3(27D 2 + 2D + 3) 



1/2 



and R > 

~ 2 



1 - 9D + J3(27D 2 + 2D + 3) 



1/2 



(19) 
(20) 



if £> < -1/3. 

In figure El we plot some curves of the surface density equation Q16|). the 
circular frequency equation (|T2|) with 6 = 0, the epicyclic frequency equation 
()13j) also with 6 = and the vertical frequency equation (|T7|) for some values 
of D and Q. The behaviour of the three first physical quantities, displayed 
as the multipole parameters are decreased, is very similar to those of the 
thick-disk model of the previous section. 
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Figure 7: Curves of a ^ = (lines without circles) and of v = (lines with 
circles) as functions of Q and R for D = 1, 0.5, 0, and -0.5. 

3 Discussion 

A Miyamoto and Nagai-like transformation was applied on the multipolar 
expansion up to the quadrupole term to produce potential-density pairs for 
flattened galaxies. These models, the first without the quadrupole term and 
the second with the quadrupole term, may be viewed as generalizations of 
the Miyamoto and Nagai models 2 and 3, respectively. The thin disk limit 
was also investigated and corresponds to generalizations of Toomre's models 
2 and 3. For each model we also calculated the velocity profile in the galactic 
plane and the epicyclic and vertical frequencies of oscillation of perturbed 
circular orbits. 

The major drawback of our models is that the density distribution is 
not a monotone decreasing function of the radial and axial coordinates for 
arbitrary values of the free parameters. Thus we imposed the condition that 
the derivatives of the density distributions with respect to R and z do not 
have critical points, except at the origin. These imposed restrictions on the 
possible ranges of the multipole moments. We found that, except in the 
thin case limit, those restrictions were also suficient to ensure non- negative 
circular velocities and also non-negative epicyclic and vertical frequencies 
(disk stability on the z = plane). 

D. Vogt thanks CAPES for financial support. P. S. Letelier thanks CNPq 
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Figure 8: (a) Curves of the surface density equation Ijl6j) . (b) the circular 
velocity equation (fT2|) with 6 = 0, (c) the epicyclic frequency equation ||T5|) 
with 6 = and the vertical frequency equation (|17j). Parameters: D = 1, 
Q = (solid lines); D = 1, Q = -0.25 (dashed lines); l) = 0.5, Q = (lines 
with points); D = 0.5, Q = —0.25 (lines with +). 
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